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Abstract 

In the framework of the product ansatz as an approximation for the two-baryon system 
we review in details the derivation of the isoscalar nucleon-nucleon spin orbit potential 
coming from the sixth order term of the extended Skyrme model. We show that the sixth 
order term contributes with a positive sign, as is the case for the Skyrme term, contrary 
to the claims of Riska and Schwesinger. Those authors considered only one part of the 
force due to the sixth order term and omitted the second part which turns out to be the 
dominant one. Our result is independent of the parameters of the model. 
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There have been several attempts ^ to extract the A^-A^ spin-orbit interaction from the 
standard Skyrme Lagrangian which includes the non-linear a model supplemented with a 
stabilizing fourth-order term in powers of the derivatives of the pion field. The calculations 
were based on the product ansatz for the two-baryon system as suggested by Skyrme 0. The 
advantage of this approximation is that, beyond its relative simplicity as compared to other two- 
baryon field configurations which can be found in the literature, it becomes exact for large A^-A^ 
separation[|. Therefore one can expect the asymptotic behaviour of the A^-A^ spin-orbit force 
derived from the Skyrme model by using the product ansatz agrees with the phenomenology. 
Unfortunately this is not the case for the isoscalar component of the spin-orbit force. All the 
authors who worked on that subject agree on the result that the standard Skyrme model predicts 
an isospin independent spin-orbit force with the wrong sign. Namely, it predicts a repulsive 
interaction while the phenomenological Bonn potential ||^ as well as the Paris potential 0] give 
an attractive one. A natural reaction that one might have in order to cure this illness, is to 
improve the Skyrme model. Indeed, it has recently been shown |0, through the study of 
some properties of the nucleon, that in order to describe properly low-energy hadron physics 
one should not restrict oneself to the standard Skyrme model but consider extensions of this 
model including higher order terms in powers of the derivatives of the pion field. Expressed in 
terms of an SU(2) matrix U which characterizes the pion field, a 6th-order term corresponding 
to a;- meson exchange [^, 



A = -^B,{u)B^{u) , (i; 



where B^' = e^'^^^Tr ( {dyU)U+{dJJ)U+{dpU)U+ ) /247r2 is the baryon current §, the uj- 
meson mass and (3^ a dimensionless parameter related to the ^ 717 width, might be a good 
candidate to solve the problem of the A^-A^ isoscalar spin-orbit force. Riska and Schwesinger [|10l 



at first and Kalbermann and Eisenberg |TT| more recently examined the influence of such term 
on the spin orbit interaction|. These authors claimed that the inclusion of the sixth-order term 
leads to the correct sign (attractive interaction) for the isoscalar spin-orbit potential. However 
they considered only one part of the interaction due to the sixth-order term in their calculations 



and omitted the second part which arises from the exchange current [10, 12 1. The aim of this 



-'^The region of validity of the product ansatz corresponds to a relative distance r larger than 1 fm 
dilaton field has also been included in llltl. 
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paper is to display the derivation of the isoscalar spin-orbit potential and show that the term 
omitted in |T0[ contributes significantly to that force. 

For a system of two interacting solitons, Skyrme |^ suggested the use of the product ansatz. 
Rotational dynamics are also introduced to obtain the appropriate spin and isospin struc- 



ture |]13|]. Thus, the field configuration of the two-nucleon system separated by a vector r can 
be written 

U{Ai,A2,^,r) = UiU2 = A,UH{ri)AtA2UH{r2)At , 
ri=x-r/2 , r2 = x + r/2 , 

where Ai and A2 are SU(2) matrices. To carry out a simultaneous quantization of the relative 
motion of the two nucleons and the rotational motion we need to treat r, Ai and A2 as collective 
coordinates. Hence we make all these parameters {r,Ai,A2) time dependent. In Eq. (H) Uh is 
the commonly used SU(2) matrix for a single soliton with the hedgehog ansatz: 

^/^(x) =exp[2T.xF(|x|) ] , x = x/|x| , (3) 

where -F(|x|) obeys the usual boundary conditions for winding number one, and the r^'s are 
the Pauli matrices. Let us observe that the chiral angle F, which will be used below in the 
numerical calculations, is obtained by solving the static Euler-Lagrange equation derived from 
the extended Skyrme Lagrangian including the non-linear a model, the 4th-order term, the 
6th-order term as well as a small chiral symmetry breaking term p|. 

The spin-orbit potential will emerge due to a coupling between the relative motion and the 
spins of the two nucleons so that we have to calculate the kinetic energy corresponding to (|l]) . 
As it is given in Ref. []TU| it reads 

K,ir) = J d'x B\U,U2) , (4) 

UJ 

where B is the spatial-component of the baryon current defined after Eq. (|I]). Before going 
further, a word of caution should be given here. In principle before identifying and extracting 
the spin-orbit potential one has to treat carefully the conversion from velocities (orbital and 
rotational) to canonical momenta ( orbital momentum and spin). Namely one has to start 
from a Lagrangian formalism, consider the "classical" kinetic energy (which is the opposite of 
Eq. (^) in the case of (|l|)), calculate the mass matrix and then invert it properly in order to 
move to a Hamiltonian formalism p[. However for a large relative distance r, the region of 



validity of the product ansatz, this procedure is equivalent to that of Refs. |]T], which 
we will use here [|T^. In the latter one starts from Eq. in our case, and make the usual 



identifications 13, 15 



-pin) i . g(n) 

"■"^ M ' ^n = --TriTAtAr.)^— , n = l,2 , (5) 

where p*^") and s*^") are respectively the radial momentum and the spin of the n-th nucleon 
while M and A are the mass and the moment of inertia of the single soliton. 

Inserting the product ansatz into the spatial components of the baryon current gives 

B([/it/2)=B(f/i)+B(f/2)+Bex(t/l,f/2) , (6) 

where B(t/„) and Bcx(t^i, U2) are the single baryon current and the exchange current densities, 
respectively. The single baryon current is given by 

j(n) g(n) \ 



B(f/„) = 5o(|rJ)(^ + r„x V) n = 1,2 , (7) 



where -Bo(^) = — -F'(r) sin^ F(r)/27r^r^ is the baryon density, while the rather complicate ex- 
pression of Bex reads 

/ 3(1) ^ pW ^(2) ^ p(2)~ 

Bi ex = DpqiC) I Aipql + A^p^i -j^ + Bipgl + B^p^i -j-j- 

in which C = AfA2 and Dpq is the 3x3 rotation matrix in the adjoint representation 

DpqiC) = ^TiiTpCrqC^) . (9) 

The sum from 1 to 3 on repeated indices in Eq. (||), and from here on, is understood. The four 
tensors appearing in the expression of Bex depend only on the positions of the two nucleons. 
They read 



(10) 



Ap.Kri, r2) = {2e.,keabp - T« a\f) /Att' , 

•^ip?/(ri, = {2tijk(^abp Rw^ Rfa R^qk " "if) /^^T^ , 

Bipql{ri,Y2) = {^—2tijkeabq T^f^ _R^p _R^f Otpi^ /4vr^ = —Aiqpl{ — Y2,~ri 

^ipqli^li^-i) = (2eijkeabq R^l Rkp ^aj " -Rgi^ <^pi) I ^'^'^ = '^iqpl{-''^2, -^l) , 
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_| I J? — T?( \ P' — — — n ^ r\ 



where 

R^i^ = -^Tr(raf/^(r„)9it//^(r„)) = r„ |^c„s„(5ia - KiKa) + y^KiKa + sleiaiKi^ , 

Tft^"^ = 2 ejjb Tni R[c^ = 2 (-s^(5bc - r„br„c) + SnCn rni) , ^^^^ 

'^m ~ ^jjfc ^abc -tl-bj ^ck ? 

9F„ _ cosF„ _ sinF„ 

By inserting now the expression in the kinetic energy we obtain 

K^{r) = I d'x (2B(f/i).B(f/2) + B^^ + 2 (B([/i) + B(t/2) ) .Be. + •••), (12) 

where we omitted the terms corresponding to the sum of the squares of the single baryon 
currents since they do not bring any couphng between the two nucleons and thus do not 
contribute to the spin-orbit force. From Eqs. and (P) one sees that the exchange current 
Box contains the isospin factor t^^Kt^"^^ due to the projection theorem 

(iV(iV^|/^a.(C)|iViiV2) = i(T«.T(2))aWaf , (13) 



while the single baryon currents are isospin independent [cf., Eq. (|^) ]. Therefore the isoscalar 
component of the spin-orbit force arises only from the first and the second term in the expres- 
sion ([T^). The calculation of the spin-orbit potential derived from the first term is straightfor- 



ward. By inserting the definition of the single currents (|^ and keeping the terms proportional 
to L.S where S = s*^^^ -|- s^^-* is the total spin and L = r x p, is the angular momentum, p being 
the relative momentum, i.e., p = p^^-* = — p*^^-*, we obtain 

/55 f.3 1 



2m,^ 

Li 

where 



/d'.2B(C/0.B(C/.) ^^_E,ML.S , (14) 



^o{r) = -/d'^ Ba{n)B„{r^) = J i'x s^sjFiFi . (15) 



The above formulae are equivalent to that obtained in Ref. [|I0[ and the function Eg as defined 
in Eq. (^) is obviously negative (see Fig.l). 

The second term in the kinetic energy (p!2| ) has been omitted in the calculations of Ref. [|10 



In fact in their paper Riska and Schwesinger ITO] referred to the article [12| for the expression 



of the exchange baryon current. However in that paper [12| the momentum dependent terms 
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as given here in Eq. (j^), have been dropped from the exchange current, since their contribution 
is of less significance than the spin dependent terms for the calculation of the deuteron form 
factors, which was the purpose of the article [|I2|; and thus no spin-orbit coupling will arise 
from their incomplete exchange current [|10] . This is obviously not the case when one takes into 
consideration the entire expression (|^). Even though the exchange current (H) is proportional 
to the rotation matrix Dpq, and thus to t^^\t^'^\ B^^ contains an isoscalar component. Indeed 
one has the formula 

Dpq{C)Dpiqi{C) = -Spp'6gq' + -ep^y eqq'q^' DpUqU^C) + C^qp' q'l^i^^n i^) , (16) 
o mn 

which can be obtained by expressing the Dpq{C) in terms of the Wigner P-functions. By 
replacing now Dpq{C)Dpiqi{C) by \5ppi5qqi (we are indeed interested only in the isoscalar part of 
the spin-orbit force) in the expression of B^^ and keeping only the relevant terms which might 
generate an L.S contribution we obtain 



2m2, 



After some tedious calculations]^ [cf. Eqs. (|l3)Q) ]; the term between brackets in the above 
equation turns out to have the structure S6cx('") ^wi so that one obtains 

^'^ ' d^x Bi ^ S6ex(r) L.S . (18) 



2m2 J 2m2 MX 

The expression of Seex reads 



/ d^x vA (if.r, (r,h)^(.? - F['){sl - F',^) + h^F^' - sl)x 



4 _ 2 2\ 
^2 S1S2) 



(2F;^ + 4 + si) + f.r, ( ^slFiF^ - \f[%- - \s\f',^ - \sIf[- - \s\f',- - \ 

(19) 

We refer to Eqs. (|3| JTT|) for the notations used in the above formula. Thus the isospin indepen- 
dent spin-orbit force generated by the 6th-order term in powers of the derivatives of the pion 
field (P reads 



^These have been checked by computer algebra. 
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FIG. 1 : The functions Eg (dashcd-dottcd line), Egox (dashed Hne) and Egtot (full line) 
in fm^^, as given in Eqs. (n3,n9l EOt), with respect to the relative N-N distance r. 



where the expressions of S6(r) and S6ex('") are given in Eqs. ([15|) and ([19|) respectively. In Fig.l 
we plot these two functions, with respect to the relative distance r. We see from that figure 
that S6(r) is negative, as it was expected and found in |10|. However the contribution of the 
exchange current which has been omitted in |T^, Sgex, is positive and larger than S6(r) so that 
the total function Sgtot, also displayed in Fig.l, is always positive. Therefore the 6th-order term 
generated by tu-meson exchange (|lD gives a repulsive isoscalar spin-orbit interaction, contrary to 
the claims of Ref. [|TU1. Our result is independent of the parameters (3^, m^^ as one can see from 



Eq. (pOl). Of course the chiral angle F on which the functions S depend, depends implicitly on 
the parameters of the extended Skyrme model [|, (/jr, e, (3uj). Nevertheless the main result of 
this paper, namely, the dominance of the positive contribution Sgex over the negative one Eg, is 
independent of the parameters of the model. We checked this result numerically by considering 
several sets of parameters which can be found in the literature; realistic as well as unrealistic 
ones. 
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In this letter we derived in details the isospin-independent spin-orbit force from the sixth- 
order order term in powers of the derivatives of the pion field (|1]). By taking into account the 
results of Refs. JI], concerning the Skyrme term [Q], we arrive to our major conclusion: neither 
the Skyrme term which can be derived from a local approximation of an effective p model, 
nor the sixth-order term generated by cu-meson exchange, reproduces the phenomenological 
isoscalar A^-A^ spin-orbit interaction. Indeed both terms give a repulsive force while it should 
be attractive according to 0, Hence the problem of the sign for this force is still unsolved 
and one should try to understand why. The key to the problem is not at hand yet but one 
can suggest some hints. The spin-orbit force is a relativistic problem and since we are in a 
region of large distance (the region of validity of the product ansatz) where this force is weak, 
any relativistic effects, even small ones, can affect that interaction. Then one should consider 
such relativistic effects in the calculations (we are thinking, e.g., to the Thomas precession 
effects due to the rotation of the nucleons). Also one sees from the Bonn potential that 
the scalar degrees of freedom, namely the a-meson which can be viewed as the responsible for 
the enhancement of the vrvr S-wave, provide one of the major contributions to the isoscalar 
spin-orbit force. So maybe one has to investigate the spin-orbit part of the two-pion exchange 



potential within the Skyrme model in order to correct the anomaly of that sign. This idea 
to account for a scalar field has been investigated in Ref. in which a dilaton field is coupled 
to Skyrmions in order to mimic the scale breaking of QCD. Even though it has been pointed 
out in Ref. [|T^ that a dilaton field is not suitable to provide a good description of low-energy 
hadron physics, the fact remains that a combination of the sixth-order term (|l|) and dilaton 
coupling, i.e., scalar meson coupling, yields an attractive isoscalar spin-orbit force as it has 



been claimed in [|Tl|. However this result remains questionable since these authors, as those of 



Ref. |]TU[, did not take into account the exchange current contribution. 
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